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Abstract: The annihilation type diagrams are difficult to calculate in any kind of models or method. Encouraged 
by the the successful calculation of pure annihilation type B decays in the perturbative QCD factorization approach, 
we calculate the pure annihilation type D — > PP{V) decays in the perturbative QCD approach based on the kr 
factorization. Although the expansion parameter l/mo is not very small, our leading order numerical results agree 
with the existing experiment data for most channels. We expect the more accurate observation from experiments, 
which can help us learn about the dynamics of D meson weak decays. 
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1 Introduction 

After decades of study, the D meson decays are still 
a hot topic in both theoretical side and experimental 
side, since they can provide useful information on flavor 
mixing, CP violation, strong interactions and even the 
new physics signal ^IHSI]- For example the recent obser- 
vation of -D" — D° mixing provides us a new platform to 
explore new physics via favor-changing neutral currents. 
By now, The CLEO-c and two B factories experiments 
have given many results about the D decays. The BES- 
III experiment is expected to give more results. The 
accurate observation can help us understand the QCD 
dynamics and the D meson weak decays. In recent years, 
many theoretical studies on the decays of D meson have 
been done based on diagrammatic approach [j] , the final- 
state interaction effects Q , combination of factoriza- 
tion and pole model 0] , factorization assisted topological 
diagrammatic approach Q, and the perturbative QCD 
(PQCD) approach 0. 

Most of the theoretical study show that the annihila- 
tion type diagrams in hadronic D decays play a very im- 
portant role jj, ■ For example in ref. [4| , the authors 
take the model-independent diagrammatic approach to 
study the two-body nonleptonic D decays, with all topo- 
logical amplitudes extracted from the experimental data. 
Their analysis indicates that the SU(3) breaking effect 
and the annihilation type contributions arc important to 
explain the experimental data. The importance of anni- 
hilation diagram contribution is also reflected from the 
large difference of and _D+ lifetime. However, these 
annihilation type diagrams are usually very difficult to 



calculate, since factorization may not work here. In 
ref. [3] , the authors use the pole model to give large anni- 
hilation diagram contributions. It is worth of mention- 
ing that the annihilation type diagrams can be perturba- 
tively calculated without parametrization in the PQCD 
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For these 



approach based on fc^ factorization 
pure annihilation type B decays, the predications in the 
PQC D appr oach have been confirmed by experiments 
later 
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The factorization that is proved in the l/rrii, expan- 
sion, can be applied to the corresponding D meson de- 
cays straightforwardly. However, the expansion is much 
poorer in D Decays than that in B decays due to smaller 
D meson mass. Anyway since there is no better method 
for the annihilation diagram calculation, the pure anni- 
hilation type decays Z?" — > K°(j) were calculated in the 
PQCD approach [9], with a good agreement with the 
experimental result. In this work, we use the PQCD ap- 
proach to analyze the 10 modes of pure annihilation type 
D —5- PP(V) decays. By keeping the intrinsic transverse 
momentum fcy of valence quarks, the end point singu- 
larity, which will spoil the perturbative calculation, can 
be regulated by Sudakov form factor and threshold re- 
summation. Therefore, the PQCD approach can give 
converging results with predictive power. 

In standard model, two body hadronic D meson weak 
decays are dominated by the contributions from tree op- 
erators, since the contributions from the penguin oper- 
ators are suppressed both by the small elements of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix and by the 
relatively small b quark mass in the c—b — u penguin 
diagram. This is in contrast to the penguin amplitude 
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in B decays, which can profit from a iarger CKM ele- 
ment and a much larger t quark mass. Although the 
suppressed penguin diagram contributions r nay be the 
main source of the direct asymmetry @, H, d, [lal , we ig- 
nore the penguin contributions in this work due to the 
small effect on the branching fractions. 

2 Formalism and Perturbative Calcula- 
tion 

For the pure annihilation type D — > PP{V) decays, 
at the quark level, the dominant contributions are de- 
scribed by the effective Hamiltonian H^ff 

H^ff = ^k,k;, [c^iti)o,{p) + c,{f,)o,{f,)] , (1) 

where V^qi and V^, are the corresponding CKM matrix 
elements, with g''' = d,s , and Ci.2(/i) are Wilson coef- 
ficients at the renormalization scale /i. Oi,2(m) are the 
four quark operators from tree diagrams 

Oi = {q'aCfi)v-A{upqa)v-A, O2 = {q'^Ca)v-A{ufiqfi)v-A, 

where a and /? are the color indices, {q'^c^)v-A = 
q'al^{^ — l^)ci3- Conventionally, the combination of Wil- 
son coefficients can be defined as 



ai=C2 + Ci/3, 02 = ^1+^2/3. 



(2) 



In the hadronic matrix element calculation, the decay 
amplitude can be factorized into soft($), hard(H), and 
harder (C) dynamics characterized by different scales 

alia, 



A ^ J dx1dx2dx-j.h1dh1b2dh2b-j.db3 

i?(x„6„t)5,(a;,)e-sW], (3) 

where hi is the conjugate space coordinate of quark's 
transverse momentum kiT, Xi is the momentum fractions 
of valence quarks, and t is the largest energy scale in the 
hard part function H{xi,hi,t). C{t) are the Wilson coef- 
ficients with resummation of the large QCD corrections 
of four quark operators. The large double logarithms 
In^Xi are summed by the threshold resummation to give 
a jet function St{xi) which smears the end-point singu- 
larities on Xi 18]. The Sudakov form factor e~'^**^ is from 
resummation of double logarithms, which suppresses the 
soft dynamics effectively and the long distance contri- 
butions in the large h region [l^ [2^. Thus it makes 
the perturbative calculation of the hard part H reliable. 
The meson wave functions $i , are nonperturbative input 
parameters but universal for all decay modes. 




M2 




Fig. 1. The diagrams contributing to the pure an- 
nihilation type D-^PP{V) decays in PQCD 



The leading order Feynman diagrams of the consid- 
ered decays are shown in Fig.l. For D — > PP decays, 
the amplitude from factorizable diagrams (a) and (b) in 
Fig.l is 

/■I /-i/A 
Aaf = -8CFfoTTm% / dx2dxz I b2db2h-idb3 
Jo Jo 

^ { [^^M2 {x2)ro2ro-i{(l)M3 (2^3) (2^3 - 2) -X30^^^ (Xg)) 

+ (l>M2i^2)<l>M^{x3){x3-l)] haf{a,^,b2,b3)Eaf{ta) 

+ [2</>M3 {x3)ro2ro3 {(fli^ {x2){x2 + {x2 ) (a;2 + 1)) 
+X2(j)M2ix2)4>M,{Xz)\ Kf{a',pM,b2)E^f{h)} , (4) 

where, Cp — 4/3 is the group factor of SU{3)c, and 
'"02(03) = "^02(03)/"^D with the chiral mass mo2(o3) of the 
pseudoscalar meson. The hard scale and the func- 
tions Eaf and h^f can be given by 



= Inax{^J{rl+X2{l-rl)){l~r'^){l~Xs)mD, 
^l-x,il-rl)mn,l/h2A/h,}, 
U = r^lSiyi{^/{l-rl){rl+X2{l-rl))'mo, 

l/&2,l/63}, (5) 

E^f{t) = a,(i)-exp[-5M,W-5M3W], (6) 



{-fHl^\pb2)Si{x,) 
e{h2^h-,)Hl^\ab2).h{ah-,) + 
e{b,-b2)Hl^\ah,)J,{ab2) \ , (7) 



with r2(3) = mM2(3)/'7iD, = (1 - 2:3(1 - r^))m|,, 
/32 = {rl + X2il - rl)){l ^ rl){l - X3)ml and a'^ = 
{rl+X2il-rl)){l~rl)ml. 

For the so called non-factorizable diagrams (c) and 
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(d) in Fig.l, the decay amplitude is 



For those D — PV decays, the decay amplitudes are 



M 



anf 



with 



16^/ -C'^Trm^ 



dxidx2dx3 



l/A 



bidbib2db2 



X 0c {Xi , 61) { [(l)tt^ (x2)Kl3 (^3) (^1 + ^2) 
+?'02'"03 {Kl2(.^2) (0M3 (2^3) (2:1+ 2^2 -2:3 +3) 

+<^M3 (2^3)(1 - 2;i + a;2 - 2^3)) 

+'/'m2(2^2) (2^3) (2:1 +2:2 +2:3-1) 

-X^-X2 + l)))] 
■hanflia,y/\M\,bl,b2)Eanfitc) 

+ro2ro3 {(j)M^ (xa) (0^/3 (2:3) (2;i - ^2 + 2:3 - 1) 

+0M3 (2;3)(2;i - 2:2 - 2:3 + 1)) 



+ 0M3(^3)(2;3- 



M3\ 

^ffllS^^) {(I)M^{X3){X2+X:i-X^' 



1) 



W2 



X^+X2~X.j,)))] 



(8) 



max{ ^J{rl+X2{l-rl)){l-rl)(l-Xs) mo , 

^/l^[{l-rl){l-X2)-x^][rl+X;{l-rl)]mD, 
l/h^Mb2), 

max{ A/|(a;i-r|-a;2(l-r|))|(l-r^)(l-a:3)mD , 
^(rl+X2{l-rl)){l-rl){l-X3)mo. 
lAi,l/&2}, (9) 



Eanf = as(t)-exp[-S'i5(t)-S'A/2(t)-5M3(t)] I 62=63' (10) 



M 



anf 



8CFfDTTm% / dx2dxz i b2db2b3db3 
Jo Jo 

X { [2(j)'^j^ {x2)ro2rv{(t)v{x3){x3 - 2) - Xscjy'yixs)) 

+<j)tt^{x2)Mx3)K~l){x3-l)] 

■haf{a,l3,b2,ba)Eaf{ta) 

- [~2(f>l.{x3)ro2rv{(t>M2 (2:2) (2:2 - 1) + (2^2 + 1) 

•0M2 (^2)) + (2:2 + (1 - 2x2yy)K, i^2)(bv {X3)] 

■Kf{a',l3,b3,b2)Eafitt)}, (12) 



[2 

IGy -^CpTTra]j J dxidx2dx3 61^6162^62 



x0o(a;i,6i) { [(j>fj^{x2)<j)v{x3) 

{xi +0:2 + (-2:1 - 2x2 + X3 + 1)?'^)) 

+ro2rv {4>1i2 (2^2) (0v (2:3) (1 - 2;i - 2:2 ^ 

+ (t>y{x3){Xi+X2~X3-l)) 

+'/'m2 (2^2 ) ((/-V (2:3 ) (2;i + 2:2 + 0:3 - 1) 
-0t'(x3)(a;i+a;2-X3 + 3)))] 

•/ian/l(Q!, VW\MM)Eauf{tc) 

-[4>i,^{x2)(Pv{x3){x3~l){2rl- 



X3) 



1) 



+n)2f'v (2^2) (0v (2^3 ) (2;i - X2 + 2:3 - 

+0V(2;3)(2;i -2;2- 2:3 + 1)) 

+'/>M2 (2^2) ((/>t-(2;3)(2;2 +2:3 - Xi - 1) 

+ (?!)V(2;3)(1 - Xl +2:2 - 2:3)))] 



1) 



■h. 



an/2 



{a,VWlKb.2)E^^f{U)], (13) 



it: 

T 



9{b,-b2)H^'^ (a6i)Jo(a62) 



+0(62-6i)i7(^' (a62)Jo(a6i) 



ITT / / 



(1) 



^o(/^6i), /3|>0, 



(11) 



where j = 1,2, (3f = 1- 
PI = jxi-rl 
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and 



[(1- ^3') (1-2:2; 
-^2(l-rg))(l 



v/(^i+2:2(l-ri))(l-ri)(l- 



l-a;i][r2 + a;3(l- 
-r|)(l -X3)rn|3 , 
X3)m]j. The ex- 



pressions of Sn[t), SM2it)j ^M^it) and iSt can be found 
in refs.[li0,[2lj. 



with Ty =r3 = my /mo- For D VP decays, the ampli- 
tudes are 



/•I /-l/A 

8CFfDT^rn% / dx2dx3 / b2db2b3db3 
Jo Jo 



X {[2(/)^(a;2)ro3rT^(0X/3 (2^3)2:3 - 
+<,^(a;3)0v(x2)((22;3-l)r^ 

■hafia,l3,b2,b3)Eaf{ta) 

- ['^Kla{^3y03rv{4'vix2)ix2-1) 

+0^ (x2 ) (x2 + 1 ) ) - 0^,3 (2;3 ) (2^2 ) (r^ ■ 
•6„/(a',/3,63,62)£'a/(4)}, 



0M3(2;3)(2;3-2)) 

-2:3+1)] 



■1)2:2] 



(14) 
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= 16 



[2 r^'^ 

V 3^^"""^^ J dxidx2dx-i J 61^6162^62 



X(/)c(a;i,6i) { [(I)f.j,^{xs)(f>v{x2) 

(xi + X2 + {-2xi - 2x2 + l)rl)) 

+ro3rv {ipij^ {x-i) ((/)^(x2)(l - ^i - 2:2-2:3) 

+(j)y{x2){l-Xi-X2'^X-i)) 

+</'m3 (2:3 ) {4>V {^2 ){Xi+X2+X:i~l) 

X2 — 2:3 

+ 3)))] 

■hanfl{0i, \/\M\lbl,h2)Eanf{tc) 
- [Kl:,{x3)(t>v{x2) 

•(l-X3 + r^(xi-a;2 + 2x3-2)) 

+ro-irv {(pli^ixz) {cj)v(x2)(l- Xi +X2-X3) 

+4>v{x2){Xl -X2-X'i + 1)) 

(2:3 ) ('/>v (2:2 ) (2:2 + 2:3 - 2:1 - 1) 

+0t'(^2)(2;i - 2:2 + 2:3 - 1)))] 

■ Knf2{a,\/\M\M,'b2)Earif{td)^ , (15) 

with Ty =r2 — my / mo ■ The form of the wave functions 
of final state pseudoscalar mesons and vector mesons can 
be found in ref.[l3l, with the different Gegenbauer mo- 
ments used in this work as 



al=0.70,af 



2tt — 1 "-4-^ 

a;f^= 0.80, af^ = 0.60, ar* 



0A5,al = 0.70,al: 



0.36, 



2^-0.10,a2^K = 0.5, 



a 



,j,--0.2,aij, = 0m,al: 



= 0.6, 



a^^=0.70,a;^.=0.6,a^^.. =0.11. 



(16) 



Since the energy release in D decays is smaller than 
that in B decays, our light meson wave functions have 
larger SU(3) breakings in D decays. For the distribu- 
tion amplitudes oi D/Ds meson, we take the same model 
as the B meson [l3| with different hadronic parameter 
a; = 0.35/0.5 for D/Dg meson. 

With the functions obtained in the above, the ampli- 
tudes of these pure annihilation decay channels can be 
given by 



+ C2M 



'^"^ anf 



.k(.)o^(.)o 



anf 



}, (17) 



G f. 



A{D' ^ if» = ^K*.K» K-4f/ + C2Xf„^,] , (18) 



A{D' ^ i^V) = %K:K4«2^ff + C2Xf,], (19) 



V2 



AiD+ ^ K+ct^) = ^K*.K. [a.A^J + C.M^J,] , (20) 



A{D,^tt+tt'') = 



'C2{M:if- 

^0, 



anf 



Gf 



v:^vma: 



■^af 



(21) 



') 



+C2{m:j, -xl; )], (22) 



T<'")-')-KC2(X::/'"' 
TMC'-^-^)]. (23) 

3 Numerical Results and Discussions 

For numerical analysis, we use the following input 
parameters: 

/^/^^ = 0.23/0. 257GeF,/A' = 0.16GeF, U ^ O-lSGeV, 
= 0.209(0. 165)Gey, f^} = 0.217(0.185)Gey, 

/^^' =0.195(0.145)Gey,/f = 0.220(0. 185)Gey, 
|T4d| =0.2252±0.00065, |Kd| = 0.9742 ±0.0002, 
iKsl = 0.97344±0.00016, |K.| = 0.2253±0.00065, 
mo^ = lAGeV.moK = l.QGeV.Ki^co = 0.375GeF. (24) 

After numerical calculation, the branching ratios of 
these decays together with experimental measurements 
(22! are listed in Table [Tl We also list the results from 
diagrammatic approach [4| and pole model 7] for com- 
parison. 

The branching ratio obtained from the analytic for- 
mulas may be sensitive to many parameters especially 
those in the meson wave function. The theoretical un- 
certainties in our calculations, shown in Table [1] are 
caused by the variation of (i) the hadronic parameters, 
such as the shape parameters and the Gegenbauer mo- 
ments in wave functions of initial and final state mesons; 
(ii) the unknown next-to-leading order QCD corrections 
and nonperturbative power corrections, characterized by 
the choice of the Kqcd = (0.375 ± 0.05) GeV and the 
variations of the factorization scales defined in eq. ([5]) and 
eq.®, respectively. 

In hadronic D decays, the SU(3) breaking effect is 
remarkable, which can be demonstrated by the decay 
channel I?" — ?> K'^K'^, with large branching ratio from 
experimental measurement. There are two kinds contri- 
butions from the quark pair dd and ss produced through 
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_D+ — tt+tt" decay, the branching ratio vanishes due to 
the exact cancelation of the contributions from uu and dd 
components. In fact, this decay is forbidden because the 
two pions can not form an s wave isospin 1 state due to 
the Bose-Einstein statistics. Any non-zero data for this 
decay may indicate the signal of new physics beyond the 
standard model. 



Table 1. Branching ratios(10 ^) for D PP{V) decays together with experimental data [S^], the recent results 
from diagrammatic approach 4] and the predictions from pole model 01 . 



decay modes 


this work 


Br (diagrammatic) 


Br(pole model) 


Br(Exp) 






,-, r,7+0.09 
'-'•^ ' -0.08 





0.3±0.1 


0.34 ±0.08 















<0.34 








8.68±0.139 


0.8±0.2 


8.34±0.65 






44+0-20 


0.29±0.22 


0.16±0.05 


<0.56 






0.541°:?° 


0.29±0.22 


0.16±0.05 


< 1.0 


D°- 




0.0121°:°°! 


0.006 ±0.005 


0.020 ±0.006 








n n2^+° °i2 

U.UZO_Q Qog 




0.020 ±0.0020 




Dt^ 


-^TV + p° 


r, 1 1 +0.87 
^■^^-0.25 




4.0±4.0 


0.2±0.12 


Dt 


— > TT + W 


n nsn+°°29 

U.UOU_o 







2.5±0.7 


Dt- 




9 1 1 +0.87 
^.^^-0.24 




4.0±4.0 





weak vertex. In SU(3) limit, the two contributions ex- 
actly cancel with each other due to the cancelation of 
the CKM matrix elements. Thus the diagrammatic ap- 
proach [3| results in zero branching ratio for this chan- 
nel. Taking the SU(3) breaking effect in account, we 
give the result in agreement with the experimental data. 
For the decay Z?° — )■ A'°0, we reproduce the result of 
ref. [3 , which agree well with the experimental data. For 



For Df — !■ 7r+p° decay, the branching ratio is larger 
than the experimental result, while for Df — > tt+cj, it is 
much smaller than the experimental result. The reason 
is that the minus sign of /9° relative to lu is compensated 
by the asymmetric space wave function of the two final 
states which are in the P-wave state. One possible so- 
lution is the soft final-state interactions as discussed in 
ref.Q. In general, the soft final state interaction should 
be important in D meson decays, because there are many 
resonance states near the D meson mass, which may 
give severe pollution to D decays calculation. We expect 
more accurate measurements from experiments such as 
LHCb and BESS-III, which can help us understand bet- 
ter the dynamics of D meson decays. 



4 Summary 

In this work, we calculate the branching ratio of the 
10 pure annihilation type -D(s) — PP{V) decays in the 
perturbative QCD factorization approach without con- 
sidering soft final states interactions. For most chan- 
nels, our results agree well with the experimental data. 
The SU(3) breaking effect is found to be remarkable, 
which can be indicated by the large branching ratio of 
D° — )• K°K° decay. We hope that the super B factories 
and BES-III can provide more accurate measurements 
for these decays, which will help us learn about the QCD 
dynamics in D meson decays and the annihilation mech- 
anism. 

We are very grateful to Yu Xin and Yu Fu-Sheng for 
helpful discussions. 
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